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The production of a pair of jets with large rapidity separation in hadron-hadron collisions, and of forward jets in deep inelastic 
scattering, can in principle be used to test the predictions of the BFKL equation. However in practice kinematic constraints lead to 
a strong suppression of BFKL effects for these processes. This is illustrated using a BFKL Monte Carlo approach.^ 

"Presented by LHO at the XXIX International Conference on High Energy Physics, Vancouver, B.C., July 23-29, 1998. 



1 Introduction 

Fixed-order QCD perturbation theory fails in some 
asymptotic regimes where large logarithms multiply the 
coupling constant. In those regimes resummation of 
the perturbation series to all orders is necessary to de- 
scribe many high-energy processes n The Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) equation El performs such a re- 
summation for virtual and real soft gluon emissions in 
dijet production at large rapidity difference in hadron- 
hadron collisions (see Figure |l|(a)) and in forward jet 
production in lepton-hadron collisions (Figure 0(b)). In 
the latter case, resummation leads to the characteris- 
tic BFKL rise in the forward jet cross section, a ~ 
(xj/xBj) x , with A = 4CU hi2a s /7r w 0.5. Similarly, in 
dijet production at hadron colliders BFKL resummation 
gives □ a subprocess cross section that increases with ra- 
pidity difference as a ~ exp(AA), where A is the rapidity 
difference of the two jets with comparable transverse mo- 
menta pti andpT2- 

Experimental studies of these processes have recently 
begun at the Tevatron pp and HERA ep colliders. Tests 
so far have been inconclusive; the data tend to lie between 
fixed-order QCD and analytic BFKL predictions. How- 
ever the applicability of these analytic BFKL solutions is 
limited by the fact that they implicitly contain integra- 
tions over arbitrary numbers of emitted gluons with ar- 
bitrarily large transverse momentum: there are no kine- 
matic constraints included. Furthermore, the implicit 
sum over emitted gluons leaves only leading-order kine- 
matics, i.e. only the momenta of the 'external' particles 
are made explicit. The absence of kinematic constraints 
and energy-momentum conservation cannot, of course, 
be reproduced in experiments. While the effects of such 
constraints are in principle sub-leading, it is desirable to 



Figure 1: Schematic representation of (a) dijet production with 
large rapidity separation A in hadron-hadron collisions, and (b) 
forward jet production in deep inelastic scattering. 



include them in predictions to be compared with experi- 
mental results. As we will see, kinematic constraints can 
affect predictions substantially. 

2 Monte Carlo Approach to BFKL Physics 

The solution to this problem of lack of kinematic con- 
straints in analytic BFKL predictions is to unfold the 
implicit sum over gluons to make the gluon sum explicit, 
and to inmlement the result in a Monte Carlo event 
generator flu. This is achieved as follows. The BFKL 
equation contains separate integrals over real and vir- 
tual emitted gluons. We can reorganize the equation 
by combining the 'unresolved' real emissions — those 
with transverse momenta below some minimum value (in 
practice chosen to be small compared to the momentum 
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threshold for measured jets) 
sions. Schematically, we have 



with the virtual emis- 



virtual 



real 



virtual ~\-real,unres. 



(i) 



reaLres. 



We perform the integration over virtual and unresolved 
real emissions analytically. The integral containing the 
resolvable real emissions is left explicit. 

We can then solve the BFKL equation by iteration, 
and we obtain a differential cross section that contains an 
explicit sum over emitted gluons along with the appropri- 
ate phase space factors. In addition, we obtain an overall 
form factor due to virtual and unresolved emissions. The 
subprocess cross section is 



da = da x ^ /„ 



(2) 



where f n is the iterated solution for n real gluons emitted 
and contains the overall form factor. It is then straight- 
forward to implement the result in a Monte Carlo event 
generator. Emitted real (resolved) gluons appear explic- 
itly, so that conservation of momentum and energy, as 
well as evaluation of parton distributions that multiply 
da, is based on exact kinematics for each event. In addi- 
tion, we iriclude the running of the strong coupling con- 
stant. SeeH for further details. 

3 Dijet Production at Hadron Colliders 

At hadron colliders, the BFKL increase in the dijet sub- 
process cross section with rapidity difference is unfor- 
tunately washed out by the falling parton distribution 
functions (pdfs). As a result, the BFKL prediction for 
the total cross section is simply a less steep falloff than 
obtained in fixed-order QCD, and tests of this prediction 
are sensitive to pdf uncertainties. A more robust pedic- 
tion is obtained by noting that the emitted gluons (cf. 
Figure |l|(a) ) give rise to a, decorrelation in azimuth be- 
tween the two leading jets&Q This decorrelation becomes 
stronger as the rapidity difference A increases and more 
gluons arc emitted. In lowest order in QCD, in contrast, 
the jets are back-to-back in azimuth and the (subprocess) 
cross section is constant, independent of A. 

This azimuthal decorrelation is illustrated in Figure ^| 
for dijet production at the Tevatron pp collider 0, with 
center of mass energy 1.8 TeV and jet transverse mo- 
mentum pt > 20 GeV. The azimuthal angle difference 
A(f) is defined such that that cos Acj) = 1 for back-to-back 
jets. The solid line shows the analytic BFK1 prediction. 
The BFKL Monte Carlo prediction is shown as crosses. 
We see that the kinematic constraints result in a less 
strong decorrelation due to suppression of emitted glu- 
ons, and we obtain improved agreement with preliminary 
measurements by the D0 collaboration shown as dia- 
monds in the figure. 



The azimuthal decorrelation can also be studied at 
the LHC pp collider 13, which has higher rapidity reach 
than the Tevatron. Figure |^ compares the decorrelation 
at the Tevatron for p? > 20 GeV (dotted curve; same as 
crosses in Fig. ||) to that at the LHC for p? > 20 GeV 
(solid curve) and px > 50 GeV (dashed curve). We 
see that at the LHC for pt > 20 GeV the decorre- 
lation is stronger and reaches to larger rapidities than 
the Tevatron. The LHC's higher center of mass energy 
(y/s = 14 TeV) relative to pr threshold allows for more 
emitted gluons, and the characteristic BFKL effects are 
more pronounced. For the perhaps more realistic LHC 
Pt threshold of 50 GeV, the kinematic suppression is 
more pronounced, but we still see a strong decorrelation. 
In all three curves we see the suppression of the decorre- 
lation by the kinematic constraints as A approaches the 
kinematic limit, where the suppression of emitted gluons 
is so strong that the curve turns over and the correlation 
begins to return. 

In addition to studying the azimuthal decorrelation, 
one can look for the BFKL rise in dijet cross section with 
rapidity difference by considering ratios of cross sections 
at different center of mass energies at fixed A. The idea 
is to cancel the pdf dependence, leaving the pure BFKL 
effect. This turns out to be rather trickyla, because the 
desired cancellations occur only at lowest order, and the 
kinematic constraints strongly affect the predicted be- 
havior, not only quantitatively but sometimes qualita- 
tively as well BO. 



4 Forward Jet Production at Lepton Hadron 
Colliders 

In deep inelastic scattering at. lepton-hadron colliders, 
the production of forward jets □ is subject to the effects 
of multiple soft gluon emission just as in dijet production 
at hadron colliders. Now the large rapidity separation 
is between the current and forward jets; see Fig. 0(b). 
The BFKL equation resums such emissions, and it is rel- 
atively straightforward to adapt the dijet formalism to 
calculate the cross section for the production of a forward 
jet with a given kx and longitudinal momentum fraction 
x.j xbj. In fact there is a direct correspondence be- 
tween the variables: px2 <-* kx and A <-> \n(x j / x^) . In 
the DIS case the variable pxi corresponds to the trans- 
verse momentum of the qq pair in the upper 'quark box' 
part of the diagram. In practice this variable is inte- 
grated with the off-shell 7*<?* — ► qq amplitude such that 
Pti ~ Q 2 ■ As a result, it is appropriate to consider val- 
ues of k T of the same order, and to consider the (formal) 
kinematic limit Xj/xBj —> oo, Q 2 fixed. In this limit we 
obtain the 'naive BFKL' prediction Oj 0t ~ (xj/xBj) x , the 
analog of ajj ~ exp(AA). 

Figure ^ shows the differential structure function 
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Figure 2: The azimuthal angle decorrelation in dijet production at the Tevatron as a function of dijet rapidity difference A, for jet transverse 
momentum pt > 20 GeV. The analytic BFKL solution is shown as a solid curve and a preliminary D0 measurement^ is shown as diamonds. 
Error bars represent statistical and uncorrelated systematic errors; correlated jet energy scale systematics are shown as an error band. 



d 2 F2/dXjdkj- as a function of XBj at HERA, with 

x 3 -=0.1, Q 2 = 50GeV 2 , Q 2 /2 < k\ < AQ 2 . (3) 

The lower dashed curve is the QCD leading-order pre- 
diction from the process j*G — > qqG, with Q — g,q, with 
no overall energy-momentum constraints. This is the 
analog of the djj — > constant prediction for dijet pro- 
duction. Note that here the parton distribution function 
at the lower end of the ladder is evaluated at x = Xj, 
independent of xb 3 - In practice, when XBj is not small 
we have x > Xj and the cross section is suppressed, as 
indicated by the lower solid curve in Fig. ^. The upper 
dashed curve is the asymptotic BFKL prediction with 
the characteristic (xj/xBj) x behavior. Finally the upper 
solid line is the prediction of the full BFKL Monte Carlo, 
including kinematic constraints and pdf dependence. We 
see a significant suppression of the cross section. Wc 
emphasise that Fig. ^ corresponds to 'illustrative' cuts 
and should not be directly compared to the experimen- 
tal measurements. Nevertheless, the BFKL-MC predic- 
tions do appear to folktw the general trend of the HI and 
ZEUS measurements!!!! A more complete study, includ- 
ing realistic experimental cuts and an assessment of the 



uncertainty in the theoretical predictions, is under way 
and will be reported elsewhereEH. 



5 Conclusions 



In summary, we have developed a BFKL Monte Carlo 
event generator that allows us to include the subleading 
effects such as kinematic constraints and running of a s . 
We have applied this Monte Carlo to dijet production at 
large rapidity separation at the Tevatron and LHC, and 
to forward jet production at HERA; the latter work is 
currently being completed. We found that kinematic con- 
straints, though nominally subleading, can be very im- 
portant. In particular they lead to suppression of gluon 
emission, which in turn suppresses some of the behavior 
that is considered to be characteristic of BFKL physics. 
It is clear therefore that reliable BFKL tests can only be 
performed using predictions that incorporate kinematic 
constraints. 
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Figure 3: The azimuthal angle decorrclation in dijet production at the Tevatron (\fs = 1.8 GeV) and LHC (y/s = 14 TeV) as a function of 
dijet rapidity difference A. Dotted curve: Tevatron, pt > 20 GeV; solid curve: LHC, Pt > 20 GeV; dashed curve: LHC, Pt > 50 GeV. 
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Figure 4: Differential structure function for forward jet production in ep collisions at HERA. The curves are described in the text. 
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